On the Effects of Dissipative Turbulence on the Narrow 
Emission-Line Ratios in Seyfert Galaxies 



S.B. Kraemer^'2, M.C. Bottorff^, & D.M. Crenshaw^ 



ABSTRACT 

We present a photoionization model study of the effects of micro-turbulence 
and dissipative heating on emission lines for number and column densities, ele- 
mental abundances, and ionizations typical for the narrow emission line regions 
(NLRs) of Seyfert galaxies. Earher studies of NLR spectra generally found good 
agreement between the observations and the model predictions for most strong 
emission lines, such as [O III] A5007, [O II] A3727, [N II] A6583, [Ne III] A3869, 
and the H and He recombination lines. Nevertheless, the strengths of lines from 
species with ionization potentials greater than that of He~''(54.4 eV), e.g. 
and Ne"*"^, were often under-predicted. Among the explanations suggested for 
these discrepancies were (selectively) enhanced elemental abundances and con- 
tributions from shock heated gas. Interestingly, the NLR lines have widths of 
several 100 km s~^, well in excess of the thermal broadening. If this is due to 
micro-turbulence, and the turbulence dissipates within the emission-line gas, the 
gas can be heated in excess of that due to photoionization. We show that the 
combined effects of turbulence and dissipative heating can strongly enhance N V 
A1240 (relative to He II A1640), while the heating alone can boost the strength 
of [Ne V] A3426. We suggest that this effect is present in the NLR, particularly 
within ~ 100 pc of the central engine. Finally, since micro-turbulence would make 
clouds robust against instabilities generated during acceleration, it is not likely 
to be a coincidence that the radially outflowing emission-line gas is turbulent. 
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1. Introduction 

Seyfert galaxies arc relatively nearby {z < 0.1), low-luminosity (bolometric luminosities 
Lboi < 10^^ erg s^^) examples of Active Galactic Nuclei. Based on the widths of their 
optical emission Unes, Seyfert galaxies are generally divided into two types (Khachikian 
& Weedman 1971). Seyfert 1 galaxies possess broad permitted lines, with full width at 
half-maxima (FWHMs) > 10^ km s~^, narrower forbidden lines, with FWHM ~ 500 km 
s~^, and optical continua which are dominated by non-stellar emission (e.g., Oke & Sargent 
1968). Seyfert 2 galaxies show only narrow emission lines and the non-stellar contribution 
to their optical continua is much weaker (Koski 1978). The broad emission lines in Seyfert 
Is vary on short timescales in response to changes in the continuum flux, which indicates 
that this emission arises in dense gas within tens of light-days of the central source (e.g. 
Peterson et al. 2004), in what is referred to as the Broad Line Region (BLR). The narrow 
lines detected in both Seyfert Is and 2s can extend ~ 1 kpc from the AGN (e.g. Schmitt 
& Kinney 1996; Schmitt et al. 2003), forming the so-called Narrow- Line Region (NLR). 
Spectropolarimetry studies revealed the presence of strongly polarized continua and broad 
permitted line emission in several Seyfert 2 galaxies. This discovery led to the unified model 
for Seyfert galaxies (Antonucci 1993), which posits that the differences between the two types 
is the result of viewing angle, with Seyfert 2 galaxies characterized by the obscuration of 
their broad-lines regions and central engines by a dusty circumnuclear torus. In this model, 
the torus collimates the ionizing radiation and the resulting illumination pattern produces a 
roughly biconical NLR (Schmitt & Kinney 1996). If the AGN is viewed at high inclination, 
the bicone axis is roughly in the plane of the sky, maximizing the projected angular size of 
the NLR. 

While historically there has been some debate about the importance of coUisional and/or 
shock ionization of the NLR (e.g., Kriss et al. 1992; Wilson & Raymond 1999), the narrow- 
ness of Radiative Recombination Continua detected in XMM-A^eu;ton spectra (Sako et al. 
2000; Kinkhabwala et al. 2002; Turner et al. 2003; Armentrout, Kraemer, & Turner 2007) is 
a clear indication that the emission-line gas is photoionized by the central source. We have 
analyzed Hubble Space Telescope (HSTj/Yami Object Spectrograph (FOS) observations of 
the NLR in the Seyfert 1 galaxies NGC 5548 (Kraemer et al. 1998a) and NGC 4395 (Kraemer 
et al. 1999), and the Seyfert 2 galaxy NGC 1068 (Kraemer, Crenshaw & Ruiz 1998b) and 
i/ST/Space Telescope Imaging Spectrograph (STIS) observations of NGC 1068 (Kraemer & 
Crenshaw 2000a, 2000b), the Seyfert 1 galaxy NGC 4151 (Nelson et al. 2000; Kraemer et 
al. 2000), and the Seyfert 2 galaxy Mrk 3 (CoUins et al. 2005). The STIS observations were 
obtained in the low- resolution/long slit mode (see Woodgate et al. 1998), with which we 
were able to spatially resolve the NLRs and explore the nature of the emission-line clouds 
as a function of distance from the central sources. Using photoionization models, we were 
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able to constrain the densities, column densities, dust/gas ratios, and ionization state of the 
emission-line gas. In general, we were able to fit nearly all of the observed emission line 
fluxes and ratios by assuming the gas was photoionized solely by the central source. Never- 
theless, we found discrepancies in the predicted strengths of a few emission lines, typically 
those from the higher ionization states (i.e., those with ionization potentials above the He II 
Lyman hmit), such as N V A124cE|, [Ne V] AA3346,3426, and [Fe VII] A6087. Usually, our 
models under-predicted the strengths of these lines by factors of ~2, although the discrep- 
ancies were often much greater for N V. The poor fit to high-ionization lines was also noted 
by Oliva (1997), and may be rectified by adjusting the elemental abundances or including 
additional ionization/excitation mechanisms (e.g. Kriss et al. 1992). In the modeling of the 
NLR emission in NGC 5548 (Kraemer et al. 1998a) and NGC 1068 (Kraemer et al. 1998b), 
we were able to improve the fit for the N V and [Ne V] lines by assuming super-solar abun- 
dances of nitrogen and neon. However, including higher abundances for selected elements 
resulted in over-predictions of lines from lower ionization states, e.g. N IV] A1486. Indeed, 
in most cases, the strengths of the lower ionization lines from these elements, e.g., [Ne HI] 
A3869 and [N II] AA 6548,6583, are consistent with roughly solar abundances (see Section 
2.1). 

In addition to photoionization model studies of the NLR, we have also used spatially 
resolved STIS spectra to analyze the kinematics of the emission-line gas in NGC 1068 (Cren- 
shaw & Kraemer 2000; Das et al. 2006), Mrk 3 (Ruiz et al. 2001, 2005) and NGC 4151 
(Crenshaw et al. 2000; Das et al. 2005). Although the kinematic studies utilized the [O III] 
A5007 line, we have found that the higher ionization lines are spatially co-located with the 
[O HI] knots (e.g. Collins et al. 2005), therefore it is probable that similar dynamical effects 
drive the higher ionization gas. Each of these sources show a similar kinematic profile, in 
which the projected radial velocity, Vr, gradually increases from the central point source 
out to r ~ 100 pc, after which the velocities abruptly begin to decrease towards systemic. 
We have been able to model the kinematics by assuming the gas is confined to a hollow, 
bi-conical envelope, whose apex is roughly coincident with the AGN. The gas is continuously 
accelerated, then begins to undergo a rapid deceleration. The dynamics of this process is 
still not clear (e.g. Everett & Murray 2006; Das, Crenshaw, & Kraemer 2007), but it is 
likely that radiation pressure from the continuum source makes an important contribution. 
The [O III] lines are also quite broad in these objects, e.g., FWHM > 1000 km s~^ near the 
apex of the bicone and > 500 km s~^ throughout the inner 200 pc of the NLR. In general. 



^For both the observations and the models, we consider the combined N V AA 1238.8, 1242.8 (N V A1240) 
and C IV AA1548.2, 1550.8 (C IV A1550), since the doublets were not resolved in the FOS and low-resolution 
STIS spectra. 
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the widths appear to decrease with distance from the AGN, although there are points of 
abrupt increases, particularly in the case of NGC 1068 (Crenshaw & Kraemer 2000; Das 
et al. 2006), which occur near the velocity turn-over points. The physical conditions that 
produce these line widths is unclear (see Section 5.). It is possible that they result from 
the super-position of multiple individual kinematic components; if so, the flow may simply 
become more chaotic at the turn-over points, broadening the profiles. However, it is also 
possible that the emission-line knots possess significant micro-turbulence. Whatever process 
slows the knots could conceivably boost the micro-turbulence. Another piece of evidence 
that the outfiowing gas is turbulent is found in the blue-shifted UV absorption detected 
along the hne-of-sight to many Seyfert 1 galaxies (Crenshaw et al. 1999; Crenshaw, Krae- 
mer, & George 2003; Dunn et al. 2007). The absorption lines nearly always have widths 
(> tens of km s~^) in excess of thermal (~ several km s"^) and, in several cases, we have 
measured FWHMs > several hundred km s^^ (e.g. Kraemer et al. 2001). The absorption 
lines appear /disappear over short timcscales (e.g. Crenshaw & Kraemer 1999; Kraemer, 
Crenshaw, & Gabel 2001; Crenshaw et al. 2003) without obvious changes in their widths, 
which is suggestive of turbulent knots passing across our line-of-sight. 

In Figures 1 - 3, we show the N V A1240/He 11 A1640, C IV A1550/He 11 A1640, and 
[Ne V] A3426/He II A1640 ratios as a function of projected radial distance, r, and FWHM. 
There is some evidence for higher values at r < 100 pc, although there are some points 
with high ionization ratios at greater radial distances for C IV and [Nc V]. More evident is 
the correlation of the line ratios with FWHM. If the broader lines are found in more highly 
excited components, as the ratios indicate, it suggests some physical connection. Therefore, 
our hypothesis is that the line widths result, at least partly, from internal micro-turbulence, 
which in turn affects the observed line ratios. 

A related issue is that the broad emission lines in AGN are smooth (e.g., Dietrich, 
Wagner, & Courvoisier 1999). If the lines are thermally broadened, there must an inordi- 
nately large number of individual clouds within the BLR. To address this problem, Bottorff 
& Ferland (2000) proposed that the clouds are internally turbulent, with micro-turbulent 
velocities that exceed their thermal velocities. Using photoionization models, they were able 
to generate smooth line profiles from a small number of clouds. Bottorff & Ferland (2002, 
hereafter BF2002) expanded on this by examining the effects of dissipative turbulence. In the 
case of non-dissipative turbulence, the turbulent motions persist, whereas, in the dissipative 
case, the motion is converted into heat. The corresponding increase in electron temperature 
will affect the emissivity of emission lines, particularly those that are coUisionally excited. 
There is evidence for dissipative heating in the diffuse interstellar medium of the Milky Way 
Galaxy (Minter & Balser 1997), hence BF2002 argued that such processes were likely to occur 
within the BLR gas in AGN. Since the thermal motions will decay due to the heating losses. 
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there must be some continual external source of turbulence. Possibilities include the intense 
radiation pressure experienced by gas close to the AGN and magnetic fields. Interestingly, 
both radiation pressure and magneto-hydrodynamic (MHD) Flows have been suggested as 
mechanisms for mass loss in AGN (see Crenshaw et al. 2003, and references therein), hence 
it is possible that forces that drive the outfiows also maintain the micro-turbulence. 

BF2002 parameterized the dissipative heating as a function of turbulence velocity, as 
follows: 

^3 

Q = rjyp^eigs cm"^ s"\ (1) 

where rjy is of order unity (see Stone, Ostriker, & Gammie 1998), Vt is the magnitude of the 
turbulence velocitjl^, p is the mass density of the gas, and D is the scale length over which 
the turbulence dissipates. BF2002 assumed that D corresponds to the physical depth of a 
BLR cloud, or ~ 10^^ cm, based on typical BLR hydrogen number densities [uh ^ 10^° 
cm~^; Davidson & Netzer 1979) and column densities (A'^^^ ~ 10^^ cm~^; Kwan & Krolik 
1981). For their models, they re-parameterized the heating rate as: the following: 

Q 2.3 X lO^^v^^ergs cm'^ s'^ (2) 

Ll3 

where is the turbulence velocity in units of 1000 km s~^, nio is the hydrogen number 
density in units of lO^'' cm~'^. Lis is the physical depth of the cloud, in units of 10^^ cm. In 
this paper, we examine the effects of dissipative turbulence in the NLR gas and have used 
equation (2) to calculate the heating. Following BF2002, we assumed that the turbulence 
dissipates over the full depth of the cloud. However, although the column densities of NLR 
tends to be much lower than the canonical BLR value (Kraemer et al. 1998a, 1998b; Kraemer 
et al. 2000, Kraemer & Crenshaw 2000a, 2000b), the physical depths are several orders of 
magnitude larger, due to the much lower densities. 

2. Photoionization Modeling 

We have calculated a set of photoionization models to investigate the effects of micro- 
turbulence and dissipative heating on the NLR emission-line ratios. The elemental abudances 
and other model input parameters are detailed in the next two subsections. 



^The turbulence velocity is equivalent to the velocity dispersion, cr(w), which is 1/2.355 times the Full 
Width Half Maximum (FWHM) of a Gaussian line profile. 
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2.1. Elemental Abundances 

The relative strengths of emission hnes strongly depend on the elemental abundances. 
Since the thermal equilibrium in the optical emission-line gas is driven largely by collisional 
cooling, the electron temperature is quite sensitive to the relative fraction of heavy atoms 
(neutral and ionized) in gas phase. For example, for heavy-element abundances more than 
several times solar values, the coUisionally excited lines can actually be weaker than for 
lower abundances due to the increased line-cooling rates. While the abundances of common 
elements such as C, O, and Ne scale with Z/Zq (where Zq indicates "solar abundances"), 
N scales as {Z/ZqY (e.g., Vila-Costas & Edmunds 1993). Hence, N/C-O-Ne can become 
quite large at super-solar abundances, which should result in the enhancement of lines such 
as N V A1240 and the N IV] A1486 multiplet compared to lines from these other elements. 

In our previous NLR studies (see Section 1), we have assumed "roughly solar" abun- 
dances (e.g. Grevesse & Anders 1989). Recently, Groves, Dopita, & Sutherland (2004) 
re-visited the issue of elemental abundances in the NLR, and incorporated the latest efforts 
to define solar abundances (e.g. Grevesse & Sauval 1998), which suggested somewhat lower 
N/H than previous studies. In order to reproduce the observed strengths of the [N II] AA 
6548,6583 lines with their photoionization models. Groves et al. had to assume twice the 
solar N/H ratio (interestingly, this is quite close to "roughly" solar value that we had used 
in our previous NLR modeling) . It is worth noting that similarly elevated abundances were 
determined for the blue-shifted UV absorbers in the Seyfert 1 galaxy Mrk 279 (Arav et al. 
2007), particularly since the UV absorbers and the emission- line gas in the inner NLR may 
be associated (Crenshaw & Kraemer 2005). Assuming all the elements have scaled together, 
this corresponds to a overall abundance of Z/Zq ^ 1.4. For this paper, we used the abun- 
dances from Asplund, Grevesse, & Sauval (2005), scaled in manner suggested by Groves et 
al. The logs of the abundances relative to H by number are: He: —1.0; C: —3.46: N: —3.92; 
O: -3.19; Ne: -4.01; Mg: -4.33; Si: -4.34; S: -4.69; Fe: -4.5. Note that the Ne/0 ratio 
is ~ 0.15, in agreement with current estimates that include the effect of solar CNNe mixing 
(see Delahaye & Pinsonneault 2006). 

Although much of the NLR gas may be dusty (Kraemer & Harrington 1986; Groves et 
al. 2004) we have not included dust in these models. The reasoning is twofold. First, there 
would be no depletion of neon and little if any of nitrogen (in the absence of ice mantles) 
onto grains. Second, we have found evidence that the dust/gas ratios in the NLR are below 
that of the ISM and, in some cases (e.g. NGC 1068: Kraemer & Crenhsaw 2000a), the gas 
in the inner NLR appears to be dust-free. 
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2.2. Model Input Parameters 

The photoionization models used for this study were generated using the Beta 5 version 
of Cloudy (Ferland et al. 1998). As per convention, the models are parameterized in terms 
of the ionization parameter, U = Q/ {Airr'^cnH) where r is the distance between the emission- 
line gas and the central source, uh is the hydrogen number density, c is the speed of light, 
and Q = j^^^y{Li,/hu) du, or the number of ionizing photons s~^ emitted by a source 
of luminosity ergs s~^ Hz~^. We assumed a plane-parallel ("slab") geometry. For the 
incident continuum, we used the spectral energy distribution which we employed in our 
modeling of NGC 4151 (Kraemer et al. 2005), which is parameterized as a broken power 
law of the form L^, oc z/" as follows: a = —1.0 for energies < 13.6 eV, a = —1.45 over the 
range 13.6 eV < hu < 0.5 keV, and a = —0.5 above 0.5 keV. We included a low energy 
cut-off at 1.24 X 10~^ eV (1 mm) and a high energy cutoff at 100 keV. While we compare our 
model predictions with emission-line spectra from other Seyfert galaxies, the results were 
not sensitive to small changes in the power-law indices or break energies, and using a single 
model SED makes comparison of the model predictions to the complete set of emission-hne 
ratios more straightforward. 

Based on our previous models, NLR gas within a few tens of pes of the AGN can be 
characterized assuming uh — 10^ cm~^. The line ratios of interest to this study do not change 
appreciably for lower values of uh, and the value used is well below the critical density for 
the level of Ne V (the upper level of the 3426 A transition; Osterbrock 1989). We 
assumed Nh — 10^^ cm~^, which is consistent with matter-bounded clouds over the range of 
ionization we used, i.e., —2 < log(f/) < 0. Based on and N^, the physical depth for the 
model clouds was 10^^ cm. The turbulence velocity, Vt, is directly command-able in Cloudy 
(Ferland et al. 1998). The dissipative heating, based on the density, physical depth and Vt 
(see Table 1), was introduced via the command-able additional heating term. We examined 
the effects of turbulence and dissipative heating for a range of Vt from (thermal broadening 
only) to 250 km s~^. 

3. Model Results 

In our models, we tested the effects of turbulence/heating on two emission lines that are 
strong and often under-predicted by photoionization models: N V A1240 and [Ne V] A3426. 
As shown in Figure 4, the He+^ zone, in which there is the strongest contribution to the 
He II recombination lines, is coincident with the and Ne"*"^ zones. Therefore, the results 
are given as line ratios relative to the He II A1640 recombination line. 
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For resonance lines such as N V A1240, there can be a contribution from photo-excitation 
to the upper level by continuum radiation. This has been suggested as the mechanism 
responsible for the strong resonance lines from He-like and H-like ions of Ne, O, and N 
detected in X-ray spectra of Seyfert galaxies (e.g. Sako et al. 2000; Kinkhabwala et al. 
2002; Armentrout et al. 2007). If the gas is turbulent, the absorption profiles arc broadened, 
which increases the number of continuum photons that the gas can absorb before becoming 
optically thick to the incident radiation, thereby increasing the contribution from photo- 
excitation. To test the effect of photo-excitation alone, we generated a set of models with 
turbulence but no dissipative heating. The results for N V are shown in Figure 5. In 
the absence of turbulence, the maximum N V/He II ratio is ~ 1.4, which occurs near an 
ionization parameter of log(f/) = —0.4. The N V/He II ratio increases with turbulence, 
reaching a value of ~ 4 for Vt = 250 km s~^. The shape of the turn-over at log(f/) > —0.4 
depends on the predicted He II A1640, which is relatively stronger at high turbulence as the 
model diverges from the Case B approximation (e.g. Osterbrock 1989). Therefore, elevated 
N V/He II ratios can be achieved for large micro-turbulence, although there is no effect on 
forbidden lines such as [Ne V] A3426. However, the range in ionization parameter for which 
N V/He II ratios > 2 are predicted is rather narrow and covers a range for which the gas 
is so highly ionized that most of the nitrogen is in higher ionization states than N"*"^. The 
N V is enhanced in this case because the smaller N"*"^ column densities result in lower N V 
optical depths. The effect of turbulence alone is not as strong at the ionization parameter 
at which the largest column density of N"*"^ is predicted, log(f/) ~ —1.3, since the line 
becomes optically thick close to the irradiated face of the slab. Also, the NLR must include 
gas at lower ionization, which still produces He II emission, hence the contribution from 
turbulent gas with high N V/He II ratios will be diluted. Therefore, photo-excitation, even 
when maximized by high turbulence, may not be sufficient to produce the highest observed 
NV /He II ratios (see Figure 1), unless the nitrogen abundances are > several times solar. 

The effect of dissipative heating scaled with micro-turbulence on the N V/He II ratio 
is shown in Figure 6. The peak ratio increases with heating, peaking at ~ 11 for vt — 250 
km s~^. Comparing these results to the models generated without heating, in addition to 
predicting high ratios over a broader range of ionization parameter, the ionization parameters 
at which N V/He II peaks drop with increased heating/turbulence. This demonstrates 
the relative contributions from photo-excitation and dissipative heating: when heating is 
included, there is a significant enhancement of the line in the more optically thick models, 
and the heating becomes the dominant effect for Vt > 200 km s"^. Including heating has a 
significant effect on the [Ne V]/He II ratio, as well, as shown in Figure 7. The slight trend to 
lower values of U for the peak ratio is due to the higher electron temperatures for the models 
with the greatest Ne+^ fraction. As these results indicate, [Ne V]/He II ratios greater than 
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unity require the heating resulting from vt > 100 km s ^. 

Since we have included dissipative heating, we examined the thermal stability of the 
gas. The thermal stability curves (e.g., Krolik, McKee, & Tarter 1981) are shown in Figure 
8. Although the unstable region, characterized by negative slope, is more pronounced for 
the models with vt — 250 km s""^, the gas will be thermally stable for the range of ionization 
parameters from which we examined emission-line ratios, i.e., \og{U) < 0.0. 

4. Comparison with Observations 

Since the under-prediction of the N V/He II and [Ne V]/He II ratios is the motivation 
for this study, we now demonstrate to what degree the addition of dissipative heating can 
reduce these discrepancies. As discussed in Section 1, kinematic studies have found that the 
[O III] A5007 lines in the NLR can be broad. For example, the brightest components in NGC 
1068 typically have FWHM of < 700 km (Das et al. 2006), while those in NGC 4151 
have FWHM of < 500 km s"^ (Das et al. 2005). Also, the largest FWHM detected among 
the strong UV absorption components in NGC 4151 is ^ 435 km s^^ (Kracmcr et al. 2001). 
In order to test the effects of micro-turbulence, we have opted to be somewhat conservative 
and have generated models with turbulence velocities of, 50, 100, and 150 km (which 
corresponds to FWHM f« 120, 235, and 350 km s~^) for comparison to the observed hne 
ratios. The values assumed for [/, Nh, and uh are as described in Section 2.2. 

In Figure 9, we compare the model predictions for N V A1240/He II A1640 versus 
[Ne V] A3426/He II A1640 against the observations. Although there are a large number of 
datapoints for which both the N V and [Ne V] hues are relatively weak, it is apparent that 
turbulence combined with dissipative heating is required for those points with N V/He II > 
1.5 and [Ne V]/Hc II > 0.5. The few points that lie outside the upper {vt = 150 km s~^) 
curve could be matched by higher turbulence/heating. In particular, note the datapoint 
for NGC 5548. In Kracmcr et al. (1998a), we had suggested that the large N V/He II and 
[Ne V]/He II ratios in the source might be the result of a high nitrogen and neon abundances. 
However, these results suggest that strong N V and [Ne V] lines are more likely indicative 
of turbulence/heating. 

In Figure 10, we show the model predictions for N V/He II versus C IV A1550/He II 
A1640. The model comparison shows that C IV/He II ratios > 4 are consistent with turbu- 
lence/heating. We also compared the N V/He II is to N IV] A1486/ He II. While the number 
of datapoints was small due to the relative weakness of N IV], combined with the effects 
of reddening, there is some indication that the N IV] may show the effects of dissipative 



-10- 



heating. 

An additional indication that the NLR gas is turbulent is that high LP /La ratios have 
been detected in some Seyferts. In Hopkins Ultraviolet Telescope spectra of NGC 1068, Kriss 
et al. (1992) measured a ratio of ~ 0.1, while, under Case B conditions (Osterbrock 1989), 
this ratio will be << 0.01. In our models, we find that for log(C/) = 1.5, without turbulence, 
Lf3/La ~ 1.4 X 10~^. For Vt — 150 km s~^, the ratio is ~ 0.01. Decreasing the column 
density by a factor of 10 increases Lf3/La by a similar factor. Hence, this strongly suggests 
that the gas in the NLR is turbulent, although the column densities of the Lyman-emitting 
gas may be somewhat lower than those assumed here, but well-within the range of values 
used in our previous NLR studies (e.g. Kraemer et al. 2000; Kraemer & Crenshaw 2000b). 

5. Discussion 

The inclusion of micro-turbulence with associated dissipative heating has the effect of 
increasing the strengths of the NV and [NeV] lines, relative to Hell. In the range of ionization 
parameter for which these lines arc strongest, i.e. —1.5 < log(f/) < —1.3, the volume heating 
rate for a turbulence velocity = 150 km s^^ is on the same order as the radiative heating rate. 
Spatially resolved spectra of the NLR of Seyferts (see Section 1) reveal that the emission- 
hne knots remain quite broad (FWHM > several hundred km s~^) throughout the inner ~ 
200 pc of the NLR. This suggests that if this is due to micro-turbulence, the turbulence is 
being continuously driven, perhaps by the contimiTim radiation from the central source. For 
example, in NGC 4151, which has a luminosity in ionizing photons of ~ 10^^ photons s~^ 
(Kraemer et al. 2005), an emission-line cloud, characterized by log(?7) = —1.5, uh = 10^ 
cm~^, and = 10^^ cm~^, would lie at r = 9.2 x 10^^ cm. If this gas formed a complete shell 
around the central source, to produce the dissipative heating rate assumed in our models 
(see Table 1) would require an energy input rate of 8.2 x 10^^ erg s~^, which is a fraction 
of the bolometric luminosity of NGC 4151, even during the low-flux state observed in 2002 
May {Lfjoi 5 x 10^^ erg s"-*^; Kraemer et al. 2005). Therefore, based on the energetics, it is 
possible that the continuum radiation can continuously drive the micro-turbulence. 

Since the emission-line gas is part of a radial outflow, one must consider possible connec- 
tions between the line widths and the means by which the gas is accelerated. Various mech- 
anisms have been suggested for the acceleration including radiation pressure (Blumenthal & 

Mathews 1975; Mathews 1976), ram pressure from winds (Weymann et al. 1982; Schiano 
1986; Smith 1993) and MHD-driven flows (Blandford & Payne 1982; Bottorff, Korista, & 
Shlosman 2000). A Blandford & Payne type MHD wind has (in cylindrical coordinates, 
r, z, 4>) a ;2-component of outflow velocity that increases monotonically, but the r- and 0- 
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velocity components can rapidly increase (e.g., in the BLR) but then decrease further out 
(e.g., in the NLR). Therefore it is possible to have large NLR line widths near the base and 
small line widths far from the base that have nothing to do with turbulence or an outwardly 
propagating disturbance. Our kinematic studies have shown that the emission-line gas is 
accelerated to velocities ~ 1000 km s^^ within the inner 100 pc of the NLR. Over the same 
distance scales, the FWHM of the lines decrease from 1000 km s~^ to a few hundred km s~^ 
(with some notable exceptions). Although this is consistent with super-position of velocity 
components, taken together with the correlation between the N V/He II, C IV/He II, and 
[Ne V]/He II ratios and FWHM of these hues, we interpret this as evidence that it is turbu- 
lence (largely) and not systematic non-radial components of the outflow that contribute to 
the hue widths. 

While a cloud of gas is being accelerated, the driving mechanism will generate insta- 
bilities within the cloud, such as Rayleigh- Taylor, for a radiatively driven cloud, or Kelvin- 
Helmholtz instabilities, for a cloud entrained in a wind. If the timescale for the instability to 
propagate through the cloud is less than the timescale for acceleration, e.g. to achieve the 
observed radial velocities, the cloud will fragment. It has been shown (e.g. Allen 1984; Smith 
1993) that fragmentation is inevitable for clouds with physical parameters similar to those 
derived from our photoionization models (which are the basis for the range in parameters 
assumed in this paper). There are only two ways to reconcile this. One way is for the clouds 
to constantly form and evaporate within a hot wind (e.g. Krohk & Kriss 2001). The other 
way is for the clouds to be robust to instability-driven fragmentation via internal turbulence 
(Allen 1984). The turbulence could be generated by multiple shocks which also accelerate 
the clouds (Marscher 1978). Interestingly, this requires turbulence velocities 0.1?;, < Vt < Vr, 
which is consistent with observations; for Vt > Vr, the clouds will evaporate. Marscher 
(1978) further notes that stable clouds must have radii > 10^^ cm, which is consistent with 
our modeling results. 

Micro-turbulence in NLR clouds will dissipate into heat on a time scale shorter than 
the NLR crossing time. This means that energy must be provided constantly to drive 
the turbulence. As mentioned above it is energetically possible for the continuum to do 
this. However transferring the energy into turbulent motions at the required rate may more 
problematic. If the outflow is an MHD wind magnetic fleld lines which thread NLR clouds 
also thread the portion of the disk from which the NLR material was originally launched. 
This provides a possible energy conduit to the NLR clouds. It is uncertain however whether 
disturbances launched from the disk could travel over a distance of 100 pc to be dissipated 
in NLR clouds. A combination of radiative driving and an MHD wind however offers the 
possibility that NLR materials can be driven against magnetic field lines thereby exciting 
disturbances which would dissipate more or less locally as heating. Such a mechanism is an 
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area for future study. At the very least the presence of magnetic fields may be necessary to 
avert overly strong shock formation in the highly turbulent gas. 

Therefore, we suggest the following scenario. The emission-line clouds form relatively 
close to the central source, perhaps in the inner few pes of the NLR. They are accelerated 
outwards by one of the mechanisms listed above. In the process, they develop significant 
micro-turbulence, which also serves to maintain them against disruption. As the clouds move 
outwards, the mechanism which generates the turbulence becomes weaker, which could be 
due the dilution of the radiation, or a drop in density of a wind. The dissipative 

heating would decrease, and there would be less enhancement of these emission-line ratios, 
as is observed (see Figures 1-3). 

6. Summciry 

We have explored the effects of micro-turbulence and associated dissipative heating on 
emission-line ratios for conditions applicable to the NLR of Seyfert galaxies. Based on our 
modeling results, we suggest the following: 

1. Micro-turbulence can increase the effect of photo-excitation, thereby boosting the 
relative strengths of resonance lines such as N V A1240. Turbulence alone will not directly 
affect collisionally excited, forbidden lines, such as [Nc V] A3426. However, we find that, if 
the turbulence dissipates over scale-lengths equal to the physical depths of the emission-line 
clouds, the resulting heating is sufficient to explain the highest [Ne V]/He II ratios detected 
in these spectra. 

2. Although we found that the elemental abundances may be slightly super-solar, when 
turbulence and heating are included, there is no need to invoke higher abundances to explain 
the observed hue ratios. 

3. Although we cannot rule out that the emission-line widths may partly be due to 
super-position of kinematic components, the correlation between the line ratios and FWHM 
suggests that the clouds arc turbulent and that the turbulence affects the relative strengths 
of the emission lines. This is further supported by the fact that including micro-turbulence of 
the same order of magnitude as the observed line widths can alleviate discrepancies between 
the models and the observations. While the turbulence appears to diminish with radial 
distance, the fact that the line widths exceed thermal throughout the inner ~ 200 pc of the 
NLR suggests that some process maintains the turbulence over large distances. 

4. Micro-turbulence can make the clouds robust to the instabihties generated during 
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cloud acceleration. This suggests that the hfetime of non-turbulent clouds in the inner NLR 
may be brief. Hence, it may not be a coincidence that the emission and absorption lines 
observed in Seyfert galaxies have widths that far exceed their thermal widths. 

In this paper, we have not attempted to fully explore the range of parameters, i.e., 
density, ionization, elemental abundances, and dust/gas ratio, that may exist in the NLR 
of Seyferts. Instead, we simply attempted to gauge the effects of turbulence/heating. In 
future work, it will be critical to examine the origins of the turbulence and develop a more 
sophisticated model for the dissipative heating, in particular regarding scale-lengths and the 
effects of internal magnetic field. In the meantime, these results suggest that it would be 
appropriate to include an approximation of turbulence/heating in photo- ionization studies 
of AGN. In a future paper, we plan to re-analyze the STIS long-slit spectra of NGC 1068 
and explore the role of turbulence in detail. 

This research has made use of NASA's Astrophysics Data System. We thank Gary Fer- 
land for his continued efforts in developing and upgrading the photoionization code Cloudy. 
We thank Bruce Woodgate and Fred Bruhweiler for useful discussions. We thank an anony- 
mous referee for their useful comments. 
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Table 1. Model Parameters" 



Turbulence Velocity (km s ^) Heating^ (ergs s ^ cm ^) 



50 2.9 X 10" 

100 2.3 X 10" 

150 7.8 X 10- 

200 1.8 X 10- 

250 3.6 X 10- 



^For the set of models, we assumed Nh — 10^^ cm ^, 
= 10^ cm-^, and the elemental abundances listed in Section 
2.1 

^For comparison, the radiative heating is 9.5 x 10-^^ ergs 
s-^ cm-^ and 1.1 x 10-^^ ergs cm-^, for \og{U) — —1.5 
and —1.0, respectively. 
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Fig. 1.— Dereddened NV A1240/He II A1640 ratios measured in HST/STIS long-slit spectra 
of NGC 1068 (Kraemer & Crenshaw 2000a, 2000b) as a function of radial distance projected 
onto the plane of the sky (top panel) and FWHM (bottom panel). The filled and unfilled 
squares are from blue- and red-shifted knots, respectively. Note the increase in the line ratios 
with increasing FWHM. 

Fig. 2.— Dereddened C IV A1550/He II A1640 ratios for NGC 1068. 

Fig. 3.— Dereddened [Ne V] A3426/He II A1640 ratios for NGC 1068. 

Fig. 4. — Selected fractional ionic abundances as a function of column density for an ioniza- 
tion parameter log(f/) = —1.0. Note that the fractions of N"*"^ and Ne+^ become negligible 
for logNn > 21.5, at which point the gas has become optically thick above the He II Lyman 
limit. 

Fig. 5. — The ratio of N V A1240/He II A1640 as a function of C/, for micro-turbulence 
velocities v = (solid line). 50 (long-dash), 100 (dash triple-dot), 150 (dash-dotted) 200 
(short-dashed), and 250 (dotted) km s^^. Dissipative heating was not included in these 
models, hence the enhancement of the N V is results solely from photo-excitation. 

Fig. 6. — Same as Figure 5, but dissipative heating has been incorporated in the models 
(see Table 1.). Note that 1) the maximum N V/He II ratio predicted is a factor of nearly 3 
times that for the turbulence-only models, and 2) the peak shifts to lower values of C/ as a 
function of v, as a result of the increased dissipative heating at high turbulence velocities. 

Fig. 7. — Same as Figure 6, showing the [Ne V] A3426/He II A1640 ratio. As noted in the 
text, for the densities and elemental abundances used in these models, [Ne V]/He II ratios 
> unity require the heating predicted for v > 100 km s"^. 

Fig. 8. — Thermal stability curves for our set of model parameters. The three curves are: 
no turbulence (sohd line), Vt = 150 km (dotted line), and Vt — 250 km (dashed line). 
The asterisks correspond to log(C/) = —1-4, at which [Ne V]/He II peaks, and the crosses 
correspond to log(C/) = 0.0 

Fig. 9.— Predicted N V A1240/He 11 A 1640 versus [Ne V] A3426/He 11 A 1640 are shown for 
models including dissipative heating for turbulence velocities 50 (dotted), 100 (dashed), and 
150 (dash-dotted) km s~^. The curves are generated for a range of ionization parameters 
—2.0 < log(C/) < 0.0, running counterclockwise, with the log(C/) values marked at several 
points for the Vt — 100 and 150 km s~^ models. The observed ratios are from our previous 
NLR studies; the symbols from the different datasets are NGC 4151, crosses (Kraemer et 
al. 2000), NGC 1068, asterisks (Kraemer & Crenshaw 2000b), Mrk 3, squares (CoUins et 
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al. 20005), NGC 5548, triangle (Kraemer et al. 1998a), and the NGC 1068 "Hot Spot", 
diamond (Kraemer & Crenshaw 2000a). 

Fig. 10.— Predicted N V A1240/He II A1640 versus C IV A1550/He II A1640 are shown for 
models including dissipative heating. Symbols are as described in Figure 9. 
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